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■ It was realized long ago that a shock propagating along the profile of decreasing density 

should accelerate [1, 2]. In compact presupernovae, like SNe Ib/c, the shock can become 
relativistic [3] and is able to produce a burst of X-ray and even 7-ray radiation [4, 5, 6]. 



The case of a peculiar type Ic SN1998bw, probably related to GRB980425, has shown that 
SNe Ib/c may have the highest explosion energy and highest production of 56 Ni among core- 
collapsing supernovae [8, 7]. Supernovae of type Ib/c are also interesting for theory due to 
problems with their light curve modeling. Their understanding may serve as a diagnostics of 
mass loss from massive stars at the latest phases of evolution. 

Numerical modeling of shock breakout in SNe Ib/c was done previously using some sim- 
' plifying approximations [9]. Our method, realized in the code Stella, allows us to get more 

IO ■ reliable predictions for the outburst. Improvements in the theory done in the current work: 
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multi-energy-group time-dependent radiation transfer, 
taking into account of (some) relativistic effects. 



A representative presupernova in our runs was a WR star built by the code kepler [10] 
(model 7A). Late light curves and spectra for this model were studied in [11]. We present 
UBV light curves found for the same model by Stella in Fig. 1. 

The Model 7A has mass 3.199 M@ (including the mass of the collapsed core). Its radius 
prior to explosion is not strictly fixed because the outer mesh zones in kepler output actually 
model the strong stellar wind and are not in hydrostatic equilibrium. So, we fixed the radius 
by hand and we got four models with radii from 0.76 up to 2 R & which are in hydrostatic 
equilibrium. Our results are summarized in the Table. 

Explosions with kepler [11] gave maximum temperature of photons T ~ 5 x 10 5 K at 
shock breakout. We have much finer mesh zoning at the edge of the star (down to lxlO -12 
Mq) and better physics and we find higher values of T. The Fig. 2 shows the difference 
between effective and color temperatures (labels 'e' and 'c', respectively, in the Table). One 
should note that the peak values of luminosity and temperature given in the table and plots 
in Fig. 2 do not contain the light travel time correction. Models, labeled as 3.5N were 
computed by D.K. Nadyozhin in 1992 with the equilibrium radiation diffusion hydrocode snv. 
Left plot in Fig. 3 shows one of his results. Dashed lines in Fig. 3 demonstrate the effect of 
averaging the light curve due the light travel time correction. 
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Figure 1: Theoretical UBV fluxes for the model 7A, exploded with E'kin = 1-47 x 10 ergs in 
comparison with observations of SN1962L of type lb. 



Parameters of shock breakouts 
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M and Ro — presupernova mass and radius, respectively, in solar units 

Ekin — the kinetic energy at infinity in 10 51 ergs 

L p and T p — the peak luminosity and temperature 

At is the width of the light curves at 1 stellar magnitude below L p 

We conclude that photon spectra at shock breakouts in SNe Ib/c peak near 3kT ~ 1 kcV 
and the X-ray fluences predicted can be observed by future missions up to tens of Mpc. 
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Figure 2: Effective and color temperatures of emerging radiation in one of the runs. 
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Figure 3: The left figure shows the results for shock breakout luminosity by the equilibrium 
radiation diffusion hydrocode SNV; the right one: a similar presupernova model, results ob- 
tained by Stella. Dashed lines show the light travel time corrected luminosity. In case of 
STELLA the correction is approximate. 



